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Abstract. I explain how to use a simple method to extract the physics of 
lattice Hamiltonian systems which are not easily analyzed by exact or other 
numerical methods. I will then use this method to establish the relation- 
ship between QCD and a special class of generalized, highly frustrated anti- 
ferromagnets. 



1 Introduction 

The title of my talk "CORE: Frustrated Magnets, Charge Fractionalization and 
QCD" might seem peculiar for a conference on light-front physics, however I hope 
to convince you that my topic is more relevant than it appears to be. Given my 
background, the examples I discuss won't be in the light-front formalism but I 
think it will be obvious that the general technique could be fruitfully adapted to 
treat the problems Simon Dalley talks about when he discusses the transverse 
lattice. 

Since COREpQ is a method for analyzing the physics of a general class of 
Hamiltonian lattice field theories it is clearly relevant to this class of problem. 
The question, "Why am I talking about frustrated magnets and charge frac- 
tionalization?" , requires a longer explanation. I must begin by telling you what 
frustrated magnets are, and then I can tell you what makes them interesting. 
Finally, I have to tell you why I am talking about charge fractionalization (why 
QCD should be obvious). 



*Presented at Light-Cone 2004, Amsterdam, 16 - 20 August 

"This work was supported by the U. S. DOE, Contract No. DE-AC03-76SF00515. 
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2 Frustrated Magnets 

The 1 + 1-dimensional Heisenberg anti-ferromagnet (HAF) is defined by Hamil- 
tonian 



The variable j in Eq^is an integer labelling the sites of a one-dimensional spatial 
lattice and I assume that — oo < j < oo. The operator s(j) = l/2er(j), where <x 
is a Pauli spin matrix, and it acts on the spin- 1/2 degree of freedom associated 
with each lattice site j. 

Intuitively, a good candidate for a mean-field approximation to the ground 
state of this system is obtained by having the average value of the spin operator 
(ip\s(j)\vp) reverse direction moving from site to site; i.e., this Hamiltonian favors 
states in which neighboring spins anti-align. 

"What is a frustrated anti-ferromagnet?". The simplest example is one where 
the frustration is geometrical in origin. For example, consider generalizing our 
one-dimensional HAF to a two-dimensional triangular lattice. Focusing on any 
one triangle we see that, if the spins on any two of the sites of the triangle anti- 
align, the third spin doesn't know what to do. Frustration can also arise when 
there are long range couplings. For example, if we add next-to-nearest neighbor 
couplings to an HAF; i.e., 



then, when a\ = «2, there is once again no way to get a low-energy pattern of 
anti-aligned spins. 

Frustrated anti-ferromagnets are interesting for several reasons. First, al- 
though they are easily defined, they are not well understood. Semi-classical com- 
putations suggest that the example defined in Eq|21 has a devil's staircase of 
phase transitions; i.e., as the ratio axja^ approaches a fixed value the system 
undergoes an infinite number of phase transitions. Another reason is that it is 
conjectured they exhibit charge fractionalization in all dimensions. This phe- 
nomenon is known to occur in 1 + 1-dimensional theories, but is not known to 
occur in system with more spatial dimensions. Basically, charge fractionalization 
means that what appears, at the level of microphysics, to be a theory of neutral 
bosons, turns out, at low energy, to be best describe in terms of a theory of 
charged fermions. 

3 CORE Basics 

As I already said, CORE is an acronym for the Contractor REnormalization 
group technique. It can be used to extract the physics of Hamiltonian lattice 
systems which are amenable neither to exact solution, nor conventional numerical 
techniques. A cartoon of the procedure followed in a CORE computation is shown 
in Fig^ where I specialize to the case of a 1-dimensional spatial lattice, where 
a spin-1/2 degree of freedom is associated to each lattice site j. My point is to 
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show how to define a procedure for truncating the original Hilbert space to an 
appropriately chosen subspace and then, how to construct a new Hamiltonian, 
defined on this subspace, which has exactly the same low-energy physics as the 
original theory. 

To construct the truncation of the Hilbert space we divide the lattice into 
finite size blocks. In the cartoon each of these blocks is assumed to contain three 
sites. Associated with site in a block there is a spin-1/2 degree of freedom so each 
block represents eight possible states. We next restrict the total Hamiltonian to 
a single block, which leads to an 8 x 8 matrix which is easily diagonalized. From 
its lowest lying eigenstates we then select a subset called the retained states, and 
define the restricted Hilbert space to be the one constructed by taking all tensor 
products of the retained states. Thus, for example, retaining the two lowest lying 
states reduces a problem on a lattice with V sites and an original 2^-dimensional 
Hilbert space to a new lattice with V/3 sites and a 2 ^^-dimensional space. 

Using the projection operator corresponding to the construction defined 
above, one computes the renormalized or effective Hamiltonian, r(H), by evalu- 
ating the following formulas. 

P{j) = \Mj))(Mj)\ P = RP(j), (3) 

3 

r(H) = lim [[T(i) 2 ]]-V2 [[T(t)HT(t)]] [[Titf]]- 1 '^ (4) 

t— >oo 

where T(t) = e~ tH and where [[O]] = POP, for any any operator O. 

Obviously, exactly evaluating this formula, which involves the exponential 
of the Hamiltonian, requires that one already knows all of the eigenvalues and 
eigenstates of the original Hamiltonian. Of course, if we know that, there is no 
point to this exercise. The CORE method is useful because there is a way to 
compute the renormalized Hamiltonian, to arbitrary accuracy, without being 
able to solve the problem exactly. This can be done because the renormalized 
Hamiltonian can always be written as a sum of finite-range connected operators; 
i.e., 

<H) = £ £ h r {B 3 ,B 3+1 , B 3+r ^). (5) 

3 r 

Here, each connected range-r operator h r (Bj, Bj+i, • • • , -Bj+ r -i) is a product of 
operators which act only on the spins in r-adjacent blocks. (A similar formula can 
be written to define the renormalized version of any other extensive operator.) 
It is remarkable that an accurate computation of the coefficient of each term in 
this expansion can be achieved by working on small sublattices. Furthermore, it 
is also true that one can get a very good approximation to the full renormalized 
Hamiltonian keeping only a few terms in the cluster expansion, say to range four 
or five. 

The first approximation to the coefficient of a range-r connected operator 
comes from evaluating Eq0] for a connected sublattice containing r-adjacent 
blocks and then subtracting from the result those pieces of t(H) already com- 
puted in the previous computations for smaller r. This coefficient is quickly 
improved by computing the higher range computations. The formulae for the 
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range-1 and range-2 connected operators are 

hiBj = r(H[Bj)) 
h 2 {Bj,B j+1 ) = T{H{B j ,B j+1 )-hi(B j )-h 1 (B j+1 ). 

The formulae for larger r are obtained by generalizing this procedure. 



(6) 
(7) 



3.1 Example 1: Heisenberg Anti-Ferromagnet 

Let's skip further generalities and see how this works for the case of the HAF 
defined in EqQ As in the cartoon in Fig^ divide the lattice into 3-site blocks. 
In this case the single block Hamiltonian is 



H{B j 



[s(3j) ■ s(3j + 1) + s(3j + 1) • s(3j + 2)] 
s(3j + 1) • [s(3j) + s(3j + 2)] 



c2 
°0- 



+1+2 



0+2 



(8) 
(9) 

(10) 



where by 5g+i+2 an d $o+2 we m ean the sum of the squares of the total spin 
operators for sites 3j,3j + 1, 3j + 2 and 3j,3j + 2 respectively. Clearly EqlTOl 
tells us that the eigenstates of the 3-site problem correspond to two spin-1/2 
multiplets and one spin-3/2 multiplet. Since the lowest lying spin-1/2 multiplet 
is the one for which 5(0 + 1 + 2) 2 = 3/4 and S(0 + 2) 2 = 2, we will use only this 
pair of degenerate states to define the space olretained states . Since these state 
have the same energy, it immediately follows that 



hi(Bj 



(11) 



To compute the range-2 contribution, h<i{B^ Bj + i), we have to solve the 6- 
site problem exactly and evaluate t(H) for the projection operator constructed 
from the four states obtained by taking tensor products of the two lowest-lying 
spin-1/2 states in each block. This computation is made quite simple if we observe 
that these four states can be recombined into one spin-0 and one spin-1 multiplet 
as follows: 

II) (12) 



TT), 



^ (I TD + +I IT)), 



1 

V2 



(ITI)-IIT)) 



(13) 



Each of these states has a definite total spin and definite total z-component of 
spin and since the total 6-site Hamiltonian commutes with these operators, it 
follows that T(t) contracts each of these states onto the unique lowest energy 
6-site state having the same quantum numbers. Thus, in the total spin basis 
t(H) will have the following diagonal form, where Eq is the energy of the lowest 
spin-0 state and E\ is the energy of the lowest lying spin-1 multiplet. 

/ E Q \ 

Ex 

Ei 

\ Ei J 
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The correct eigenvalues are 

E- = -2.493577... ; E x = -2.001995 .. . . (14) 

Following the approach we just used it is clear that, in the original basis of tensor 
products of retained states, this Hamiltonian takes the form 

r(H(Bj, B j+1 ) = [« 1 + ■ *ti + 1)] - (15) 

3 

where a = (E + 3E 1 )/4 = -3.99507 and f3 = E x - E = .491582. 

Since, after one step the new effective Hamiltonian is a multiple of the unit 
matrix plus [3 times the original Hamiltonian, we see that this HAF Hamiltonian, 
in range-2 approximation, is at a fixed point of this renormalization group trans- 
formation. In other words, no matter how many times we repeat this process we 
always get a Hamiltonian of the form 

r(H) n = ^Kl + /f s(j) ■ s(j + 1)] . (16) 

3 

Moreover, since (3 < 1, we see that the interaction term in the Hamiltonian 
eventually iterates to zero, which tells us that we are dealing with a massless 
theory. 

The observation that [3 n — ► tells us that all we have to do is compute the 
limiting value of a n divided by the number of sites on the initial lattice to obtain 
the ground-state energy density. Now, since after the first renormalization group 
transformation the lattice has V/3 sites, the total effect of the term proportional 
to a is to contribute a V/3 to the energy of all states in the new effective theory. 
Thus, dividing by V we obtain a contribution of a/3 to the ground state energy 
density. At this point the simplest thing to do is subtract the constant term 
from the Hamiltonian and do another renormalization group transformation. 
This time the new constant term will be /3 a and will correspond to a theory 
defined on a lattice with V/3 2 sites and so its contribution to the energy density 
will be Pa/3 2 . Proceeding in this manner we see that the total energy density 
is given by the geometric series 

f = I?(f) = 3 (1-/3/3)- 

Plugging in the values of a and (3 obtained from the 6-site computation we arrive 
at the results 

£rcn-grou P = -0.448446... ; £ cxact = -0.442147... (18) 



which shows that the simple range-2 approximation is accurate to about 1%. Not 
bad for a simple first principles computation based upon diagonalizing at most 
4x4 matrices. 
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3.2 Example 2: The 1 + 1-dimensional Ising Model 

Now that I have shown you that even the crudest approximation to the exact 
CORE transformation gives qualitatively good results for the Heisenberg anti- 
ferromagnet I wish to spend a few moments showing you the kind of accuracy 
one can obtain by working a little harder. To do this I consider the case of the 
1 + 1-dimensional Ising model in a transverse field; i.e., the theory defined by the 
Hamiltonian 



This theory is interesting for several reasons. First, it clearly undergoes a 
quantum phase transition as A varies from A = to A = ir/2. To see this observe 
that for A = the theory has a unique ground state; namely, the one which is 
a product of eigenstates of cr 2 (j) corresponding to the eigenvalue +1. However, 
for A = 7r/2 there are two possible ground states made by taking the product of 
eigenstates of <J x {j). The first, is a product of eigenstates where the eigenvalue 
of a x (j) are an equal to +1; the second, where all the eigenvalues are equal to 
— 1. This tells us that as A varies the system goes from having a unique ground 
state, to having a doubly degenerate ground state corresponding to spontaneous 
breaking of the symmetry cr x (j) —&x{j)- Even more interesting is the fact 
that the low-lying excitations of the theory for small A correspond to having 
one spin flipped from up to down, but the excitations in the case of A near 
7r/2 are kinks or anti- kinks (i.e., states where half of the system seems to be in 
one ground state while the rest is in the other). I will now show that a simple 
CORE computation is not only able to compute the ground state energy density 
to high accuracy, but it correctly finds the nature of the excitation spectrum 
and computes both the mass gap and magnetization (i.e., the mean value of the 
operator (1/V) J2{cr x (j)) as a function of sin(A). 

Figs l2l4l plot results obtained by carrying out range-3 CORE computations 
for specific values of sin(A). Range-3 computations begin by dividing the lattice 
into 3-site blocks, and retaining the two lowest lying eigenstates of the 3-site 
Hamiltonian and then solving the 6-site and 9-site problems exactly. Given the 
accuracy of the results it is remarkable that the hardest computation one does is 
to diagonalize a 512 x 512 matrix. FigElplots the difference between the exact and 
CORE values for the ground state energy density divided by the exact energy 
density for a range of values of sin(A). This is done to make the errors visible. The 
dotted curve corresponds to results obtained for a range-2 (6-site) computation, 
whereas the solid curve gives the results for a range-3 computation. Clearly 
the method makes its biggest errors near the phase transition, sin(A) = .5, but 
even there it does well. In Fig^l we see a plot of the exact mass gap (the solid 
curve) against the CORE approximation to the mass gap (open circles). The only 
significant errors on this curve are near the critical point and they are essentially 
due to the fact that at range-3 CORE makes a small mistake in the location of 
the critical point. Fig|5]is a plot of the exact magnetization (solid curve) above 
the phase transition versus the values computed by CORE on a A-by-A basis. It 
is notoriously difficult to compute this curve in a simple way since the behavior 




(19) 
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corresponds to a critical exponent of 1/8. The final figure, Fig® shows how one 
can extract both the critical point and the exponent for the magnetization by 
varying both to produce a straight line. 



4 Mapping One Theory To Another 

At his point I have shown how to use CORE to extract the physics of a lattice 
Hamiltonian theory to high accuracy. In the remainder of this talk I will focus 
on what one can learn by using CORE to map a theory into a very different 
looking but equivalent Hamiltonian system. 



4.1 Example: Massless Bose Free Field 



Let us first consider the case of a massless Bose free field in 1 + 1-dimensions 
whose Hamiltonian is 



~p(j? + ^-x(jf + (x(j + i)- x (j)y 



(20) 



My objective is to show that this theory can be mapped into an equivalent spin 
system by truncating the single-site theory to the two lowest lying oscillator 
states. Furthermore, I will show that one can obtain a very good approximation 
to the exact renormalized Hamiltonian by including only a few terms in the 
cluster expansion. 

I begin by defining my blocks to consist of a single lattice site, in which case 
the one-block Hamiltonian becomes 



En 



1 



2V^ 2 + 2 



(21) 
(22) 



where E$ is the ground state energy of this system. Since this is just the Hamilto- 
nian for a simple harmonic oscillator, setting /i = 0, we find the range-1 connected 
term is 

^ = E[-^ + ^c?)1- ( 23 ) 



Carrying out the 2-site and 3-site CORE computations we get the following 
range-2 and range-3 connected operators 



h 2 = Y J [1.1351^ + .278<r 2 (j) - mSo x (j)a x (j + 1) 
j 

-.268 cT y (j)a y (j + 1) - .183 a z (j)a z (j + 1)] 



(24) 



and 



h 3 = J2 [1-0451,- + .138a z (j) - .117 a x (j) a x {j + 1) - .202a y (j) a y (j + 1) 
j 

-.167(7,0') °*U + 1) + ■003o- a .(j) a x (j + 1) a z (j + 2) 
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-.019 ^.(j) a z (j + 1) a x {j + 2) + .003a z (j) a x (j + 1) a x (j + 2) 
+.032^0") a y (j + 1) a z (j + 2) - .019 a v (J) a z (j + 1) a y (j + 2) 
+.032 a z {j) a y (j + l)a y {j + 2) (25) 

Clearly the typical coefficients of range-2 operators are significantly larger than 
the coefficients of range-3 operators. Nevertheless, since as the range of the con- 
nected operators grow more operators appear, it is important to be sure that the 
effects of the many small terms don't overwhelm the larger terms in the renor- 
malized Hamiltonian. I show that this is not a problem in Figs EE To create 
these plots I computed the renormalized Hamiltonian out to range-7. The dif- 
ferent colored plots exhibit the result of exactly diagonalizing the Hamiltonian 
defined by keeping all terms up to range-n. (What is plotted are the exact eigen- 
values in increasing order versus the integer which labels the eigenvalue after 
sorting them in increasing order.) For this reason the last curve in FigEl and 
Fig0 represents the exact eigenvalues of the lowest lying 6-site (or 7-site) states 
which have an overlap with the retained states. As you can see, after range-4, the 
longer range terms play only a small role in determining the eigenvalues over the 
entire range of energies. In figures FigEland FigEJwe see the same plots for up 
to 10-site sublattices. Once again we see that the convergence is rapid, although 
I don't bother to plot the corresponding exact eigenvalues for 8 or 10-sites since 
it would be hard to see the differences. 

The lesson to be learned from these plots is that even in the case of what 
would seem to be a very bad approximation to a massless theory, where the 
correlation length is infinite, the finite range cluster approximation to the CORE 
transformation is rapidly convergent. 



4.2 Example: Massless Free Fermion 

Now that we have seen that a massive and massless free boson theory can be 
mapped into a spin system, with no loss in information about the low-energy 
theory, I want to do the same thing for the case of a massless free fermion. In 
this case, the Hamiltonian is 

31 ,32 

s'(ji-h) = S'{j2-h)- (27) 

In what follows I will show that fermionic theory maps onto a highly frustrated 
anti-ferromagnet, where the fundamental spin on each site correspond to single- 
site states of the original theory having zero-charge [2]. To do this I expand tpU) 
in terms of single-site annihilation and creation operators and define the four 
possible single-site states to be: 

|0> |+) I-) |±>- (28) 

Here, |0), is the state annihilated by the particle and anti-particle destruction 
operators. The states |+) and |— ) are one particle and one anti-particle states 
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(having positive and negative charge), and the state |±) is the zero-charge state 
containing both a particle and an anti-particle. Since the Hamiltonian H, Eg 1271 
contains no terms which refer to a single site, it follows that at the single-site 
level these four states are degenerate; therfore 



h = 0. 



(29) 



Since there is no reason based on energy considerations to select any two states 
to keep, we will define the space of retained states to be those generated by 
taking tensor products of the chargeless single-site states |0) and |±). With this 
choice we carry out a CORE transformation to obtain the following renormalized 
Hamiltonian 



h 2 
h 3 



E 



i 



lj + ■ *(j + 1) 



(30) 



H r , 



[-.28033!j + .9428er(j) • a(j + 1) + .2357 cr{j) ■ a(j + 2)] (31) 

3 

[-.31099!j + .80001, tr(j) ■ a(j + 1) + .23492 a (J) ■ a(j + 2) 
-.01915 a(j) ■ a(j + 3) + .03559 tr(j) ■ <x(j + 1) a(j + 2) • a(j + 3) 
-.08033 • a(j + 2) a(j + 1) • a(j + 3) 

+.0403 a(j) ■ <r(j + 3) a(j + 1) • <r(j + 2)] (32) 
E hr (33) 



As advertised, we see that the structure of the effective Hamiltonian is indeed 
that of a generalized frustrated anti-ferromagnet and, once again, the cluster 
expansion is rapidly convergent. Moreover, although this computation was done 
for the case of a 1 + 1 dimensional theory, out to range-3 the same general struc- 
ture would be obtained for a 1 + 3-dimensional theory. This strongly implies, at 
least for special values of the couplings, that a highly frustrated anti-ferromagnet 
can exhibit charge-fractionalization in higher dimensions. I say this because the 
underlying degrees of freedom of this theory are neutral, however the theory is 
equivalent to a theory of free fermions. What other surprises can lurk in this 
class of theories? 



4.3 The Lattice Schwinger Model 

The lattice Schwinger model is a 1 + 1-dimensional gauge theory defined by the 
Hamiltonian 



H = Hf — ^2 Pn \ n - m \ Pm, 

m,n 

where the charge density operator p n is defined to be 

Pn = 



(34) 



(35) 



This model is interesting in that it provides an example of a confining gauge 
theory. Moreover, the strong coupling (i.e., large e) properties of this theory are 
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very similar to the properties of strong coupling lattice QCD and, in the limit of 
large e, this problem is very easy to analyze. Finally, from the point of view of 
CORE, although this problem is very difficult to analyze using earlier real-space 
renormalization group methods there is no problem defining CORE transforma- 
tions starting from locally gauge-invariant states. This is because, as we saw in 
the free fermion theory, it is possible to restrict to locally chargeless states and 
still get a non-trivial renormalized Hamiltonian. In fact, there is even more mo- 
tivation to truncate to the single-site chargeless states, |0) and |±), which are 
the ultimate in locally gauge invariant states, since in this model they span the 
set of all low lying physical states for large e. 

The space of states spanned by the tensor products of these chargeless states 
becomes completely degenerate in the limit e — * oo because they contain no flux. 
Thus, it is possible to expand the Hamiltonian of this system about the large e 
limit by doing degenerate perturbation theory in the kinetic term. This analysis 
was carried out in Ref.[3|, where it was shown that this strong coupling theory 
was equivalent to a frustrated anti-ferromagnet. The question which dominated 
this analysis was, "To what degree would the strong coupling discussion work as 
one moves towards weak coupling?". Kirill Melnikov and I@j recently analyzed 
the lattice Schwinger model and showed that the strong coupling theory connects 
smoothly to the weak coupling theory and one can easily understand all of the 
properties of the continuum limit using a generalized form of SLAC fermionic 
derivative. Space doesn't permit a full discussion of this issue now, but I mention 
it to show that building a CORE truncation procedure using these locally gauge 
invariant states is reasonable. 

Obviously the strong coupling expansion and CORE computation agree, how- 
ever, as one moves away from strong coupling the coefficients of the range-r terms 
become non-trivial functions of e. It is important to point out, however, that the 
general structure of the range-2, range-3 and range-4 terms in the perturbation 
expansion of the strong coupling theory and corresponding CORE computation 
will, for symmetry reasons, be the same. 

Having said this we see that in the strong coupling limit the Schwinger model 
is certainly equivalent to a frustrated anti-ferromagnet and moreover, the CORE 
computation to range-3 says this will be true for all couplings. Conversely, we 
see that, at least for special couplings, the general frustrated anti- ferromagnetic 
systems can be equivalent either to free massless fermions or to fermions in- 
teracting with gauge-fields. So now we see that a frustrated anti-ferromagnet 
can be equivalent to a theory of free fermions, or fermions interacting through 
long-range gauge-fields, at least for some values of the couplings. 

I would now like to finish by talking about QCD. 

5 What About QCD? 

Lattice QCD is much like the lattice Schwinger model in that it is a theory of 
quarks interacting with color gauge-fields (defined on lattice links). Moreover, 
due to the non-abelian nature of the gauge-fields, the theory confines at strong 
coupling, in the same way fermions are confined in the Schwinger model. If one 
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uses any derivative which allows quarks and anti-quarks to live on the same sites 
then, as in the Schwinger model, one immediately finds that at infinite coupling 
the ground state of the system is highly degenerate, since all single site color 
singlet states have vanishing energy. In other words, at strong coupling quark 
states with the quantum numbers of baryons or mesons all have zero-energy. This 
means that if one chooses the space of retained states to be the one spanned by 
tensor products of these zero-energy single-site states, one can then use CORE 
to construct an effective theory of interacting mesons and baryons. If, as in 
the Schwinger model, one computes the effective strong coupling theoryjSj to 
order 1/g 2 (where g is the gauge coupling constant) one obtains a generalized, 
frustrated anti-ferromagnet whose Hamiltonian takes the form 

ffei = iE Q aa W ■ + A») SM- (36) 

For n flavors of quarks, the operators Q a (j) are the generators of the group 
SU(4n). In particular, for the case of 3 flavors the nearest-neighbor anti- 
ferromagnetic interaction leads to a theory with an SU(6) x SU(6) chiral sym- 
metry, which breaks spontaneously to produce a large number of would-be Gold- 
stones bosons. The next-to- nearest neighbor frustration terms break this sym- 
metry to chiral SU(3) x SU(3), which breaks spontaneously to a theory in which 
the vector SU(3) symmetry is exact and the axial-vector symmetry is realized 
by the existence of eight Goldstone bosons. 

The explicit form of the next-to-nearest neighbor terms in the effective Hamil- 
tonian, Eq 1361 allows us to compute the transformation properties of the SU (12) 
symmetry breaking terms. Therefore, it is possible to predict the pattern of mass- 
splittings in the large multiplet of would-be Goldstone bosons without solving 
the theory exactly. The table in Fig llUI shows the quantum numbers of these par- 
ticles and the particles with increasing mass are those shown in darker shades 
of blue. Of course one has to do more than study the Hamiltonian obtained 
from strong coupling perturbation theory to correctly understand the details of 
the masses, magnetic moments, etc.. Nevertheless, the observation that out to 
range-3 the effective Hamiltonian obtained from a CORE computation and that 
obtained from strong coupling perturbation theory will be the same, says that 
the most important terms in the effective Hamiltonian will have the same sym- 
metry structure. The difference will be that the values of the coupling constants 
appearing in front of the different terms will now be non-trivial functions of 
the gauge-coupling g. Assuming, that as in the Schwinger model, the pattern of 
symmetry breaking is the same as it is at strong coupling, one obtains several 
interesting phenomenological predictions. 

First, one obtains the usual Gell-Mann Okubo sum rules for the mass- 
splittings in the various SU(3) multiplets. Second, one obtains the additional 
prediction 

m l* -m 2 p = m 2 L -ml (37) 
which is good to about 16%, and the prediction 



mj « 1.2 ± 0.2GeV 



(38) 
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which is a bit low to match a known state, but after all, this is just a first order 
perturbation theory computation. 

In addition to these results, the approximate SU(6) x SU (6) symmetry leads 
to the famous prediction for the ratio of the proton to the neutron magnetic 
moment 

!± = -f. (39) 

2 

Even more interesting is the fact that, since the weak axial- vector current is 
not a generator of the SU(6) x SU (6) symmetry, the infamous incorrect SU(6) 
prediction that 

«± = J (40) 
9V 3 



is replaced by 



9 -± = -U (41) 
9V 3 



where X is some reduced matrix element. Combining this with the fact that the 
usual chiral SU (3) x SU (3) symmetry is exact, r the Adler-Weissberger prediction 
for this ratio should be true. 



6 Conclusions 

To summarize, I hope that I convinced you that CORE techniques can be ef- 
fectively used to study Hamiltonian lattice systems. In addition, I hope to have 
convinced you that they can be used to take lattice QCD and map it into an 
effective theory whose degrees of freedom only have the quantum numbers of 
ordinary mesons and baryons (and glue-balls, if one chooses the fundamental 
block to be a square). If one follows Simon Dalley's approach, the usefulness of 
CORE would be that this effective theory would be constructed using DLCQ 
to solve the 1 + 1-theory to high accuracy and construct the space of retained 
states. Having done this one would then couple pairs of lines, etc. and try to 
evaluate the CORE formula using similar DLCQ techniques for the coupled sys- 
tem. In this way smallish DLCQ computations could be used to construct the 
an effective Hamiltonian which contains the low energy physics of QCD. 

In any event, even if this approach doesn't appeal to you, I would suggest 
that we have at least learned something; namely, that the physics of generalized 
frustrated systems is very rich. This is because, by using CORE to map known 
theories into such systems for special values of the coupling constants, we see 
that these theories must exhibit a plethora of interesting phases. 
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Figure 1. Subdivision of 1-dimensional lattice and truncation scheme. 
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Figure 2. The difference between the exact ground state 
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Figure 3. The exact mass-gap (solid curve) versus point- 
by-point CORE computation (open circles) plotted as a 
function of A 
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Figure 4. Exact magnetization versus (solid curve) versus Figure 5. Extraction of the critical point and critical 
point-by-point CORE computation (open circles) plottedas exponent for the magnetization by fitting to the indicated 
a function of A form 




Figure 6. 6-Site Comparison 



Figure 7. 7-Site Comparison 



16 



LC 2004 




Figure 8. 8-Site Comparison Figure 9. 10-Site Comparison 
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Figure 10. Exact and approximate Goldstone bosons when the effects of next-to-nearest 
neighbor terms are included. 



